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Abstract 

Cognitive radio that supports a secondary and opportunistic access to licensed spectrum shows great 
potential to dramatically improve spectrum utilization. Spectrum sensing performed by secondary users to 
detect unoccupied spectrum bands, is a key enabling technique for cognitive radio. This paper proposes a 
truncated sequential spectrum sensing scheme, namely the sequential shifted chi-square test (SSCT). The 
SSCT has a simple test statistic and does not rely on any deterministic knowledge about primary signals. 
As figures of merit, the exact false-alarm probability is derived, and the miss-detection probability as 
well as the average sample number (ASN) are evaluated by using a numerical integration algorithm. 
Corroborating numerical examples show that, in comparison with fixed-sample size detection schemes 
such as energy detection, the SSCT delivers considerable reduction on the ASN while maintaining a 
comparable detection performance. 
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I. Introduction 

Most radio frequency spectrum is allocated primarily based on fixed spectrum allocation strategies that 
grant licensed users to exclusively use specific frequency bands to avoid interference. Recent reports [1] 
released by the Federal Communications Commission show that, a large amount of allocated spectrum 
particularly television bands, is substantially under-utilized most of the time whereas a small portion 
of spectrum bands such as cellular bands, experience increasingly congestion and scarcity due to rapid 
deployment of various wireless services. Cognitive radio, which enables secondary (unlicensed) users to 
access licensed spectrum bands not being currently occupied, can fundamentally alter this unbalanced 
spectrum usage and therefore can dramatically improve spectrum utilization. Since licensed (primary) 
users are prior to unlicensed (secondary) users in utilizing spectrum, the secondary and opportunistic 
access to licensed spectrum bands is only allowed to have negligible probability of deteriorating the 
quality of service (QoS) of primary users (PUs). Spectrum sensing performed by secondary users (SUs) 
to detect the unoccupied frequency bands, is the key enabling technique to meet this requirement, thereby 
receiving considerable amount of research interest recently [2]-[5]. 

Albeit in essence a conventional signal detection problem, the design of a spectrum sensing scheme 
needs to cope with several critical challenges that stem from special attributes of cognitive radio networks. 
First, it is often difficult for SUs in a cognitive radio network to acquire complete or even partial 
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knowledge about primary signals. Secondly, SUs need to be able to quickly detect primary signals 
at a fairly low detection signal-to-noise ratio (SNR) level with low detection error probabilities. To this 
end, several spectrum sensing schemes, such as matched-filter detection [6], energy detection [2], [7], 
and cyclostationary detection [8], have been proposed and investigated. Among these sensing schemes, 
energy detection is particularly appealing as it does not rely on any deterministic knowledge of the 
primary signals and has low implementation complexity. However, when the detection SNR is low, energy 
detection entails a large amount of sensing time to ensure high detection accuracy, e.g., the sensing time is 
inversely proportional to the square of SNR [9]. To overcome this shortcoming, several sensing schemes 
based on the sequential probability ratio test (SPRT) have been proposed under various cognitive radio 
settings [4], [10], [11]. 

The SPRT has been widely used in many scientific and engineering fields since it was introduced by 
Wald [12] in 1940s . Perhaps, the most remarkable character of the SPRT is that, for given detection error 
probabilities, the SPRT requires the smallest average sample number (ASN) for testing simple hypotheses 
[13]. In comparison to fixed-sample-size sensing schemes, the sensing scheme based on the SPRT requires 
much reduced sensing time on the average while maintaining the same detection performance [10]. 
Nonetheless, the existing SPRT based sensing scheme [10] suffers from several potential drawbacks. 
First, when primary signals are taken from a finite alphabet, the test statistic involves a special function, 
which incurs high implementation complexity [4], [10]. Second, evaluating the probability ratio requires 
deterministic knowledge or statistical distribution of certain parameters of the primary signals. Acquiring 
such deterministic information or statistical distribution is practically difficult in general. Thirdly, the 
existing SPRT based sensing scheme adopts the Wald's choice on the thresholds [12]. However, the 
Wald's choice, which works well for the non-truncated SPRT, increases error probabilities when applied 
to the truncated SPRT. 

In this paper, we propose a truncated sequential spectrum sensing scheme, namely the sequential 
shifted chi-square test (SSCT). The SSCT possesses several attractive features: 1) Like energy detection, 
the SSCT only requires the knowledge on noise power and does not rely on any deterministic knowledge 
about primary signals; 2) compared to fixed-sample-size detection such as energy detection, the SSCT is 
capable of delivering considerable reduction on the average sensing time while maintaining a comparable 
detection performance; 3) in comparison with the SPRT based sensing scheme [10], the SSCT has a much 
simpler test statistic and thus has lower implementation complexity; 4) and the SSCT offers desirable 
flexibility to strike a trade-off between detection performance and sensing time when the operating SNR 
is higher than the minimum detection SNR. To evaluate the detection performance of the SSCT, we derive 
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the exact false-alarm probability, and employ a numerical integration algorithm from [14] to compute the 
miss-detection probability and the ASN in a recursive manner. Notably, the problem of evaluating the 
false-alarm probability of the SSCT resembles the exact operating characteristic (OC) evaluation problem 
associated with truncated sequential life tests involving the exponential distribution [15] [16]. The latter 
problem has been solved by Woodall and Kurkjian [15]. Despite the similarity of these two problems, 
the Woodall-Kurkjian approach cannot be applied to directly evaluate the false-alarm probability of the 
SSCT. In addition, the Woodall-Kurkjian approach is not applicable to evaluate the ASN for truncated 
sequential life tests involving exponential distribution [16]. As a byproduct, our approach to evaluating the 
false-alarm probability of the SSCT, can be readily modified to evaluate the ASN for truncated sequential 
life tests in the exponential case. 

The remainder of this paper is organized as follows. Section [TT] presents the problem formulation 
and provides necessary preliminaries on energy detection and a SPRT based sensing scheme. Section 
UlTI introduces the SSCT and its equivalent test procedure. Section [TV] deals with the evaluation of the 
error probabilities of the SSCT. In particular, this section provides an exact result for the false-alarm 
probability, and a numerical integration algorithm to recursively compute the miss-detection probability. 
Section [V] presents an evaluation result on the ASN of the SSCT, while Section [VTJ provides several 
numerical examples. Finally, Section IVIII concludes the paper. 

The following notation is used in this paper. Boldface upper and lower case letters are used to denote 
matrices and vectors, respectively; Ij. denotes a k x k identity matrix; E[-] denotes the expectation 
operator. (-) T denotes the transpose operation; 1*. denotes afcxl vector whose entries are all ones; Hp 
denotes a set of consecutive integers from p to q, i.e., Hp := {p,p + 1, . . . , q}, where p is a non-negative 
integer and q is a positive integer or infinity; (-) c denotes a complement of a set; I{ x >t} denotes an 
indicator function defined as h x > t \ = 1 if x > t and I{ x > t } = if x < t, where x is a variable and t is 
a constant. 

II. Problem Formulation and Preliminaries 

In this section, we start by presenting a statistical formulation of the spectrum sensing problem for 
a single SU cognitive radio system. We next give a brief overview on two sensing schemes, energy 
detection and a SPRT based sensing scheme, which are closely related to the SSCT. 
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A. Problem Formulation 

Consider a narrow-band cognitive radio communication system having a single SU. The SU shares 
the same spectrum with a single PU and needs to detect the presence/absence of the PU. Let Ho and 
Hi denote the null and alternative hypotheses, respectively. The detection of the primary signals can be 
formulated as a binary hypothesis testing problem as follows 

H :ri = Wi, i = 1,2,... , (1) 
H 1 :r i = hsi +w i: i = 1, 2, . . . , (2) 

where is the received signal at the SU, wi is additive white Gaussian noise, h is the channel gain 

between the PU and the SU, and Si is the transmitted signal of the PU. We further assume that 

Al) Wi's are modeled as independent and identically distributed (i.i.d.) zero mean complex Gaussian 

random variables (RVs) with variance <r^/2 per dimension, i.e., Wi ~ CM(0, <t^), 
A2) the channel gain h is constant during the sensing period, 
A3) the primary signal samples Si are i.i.d., 
A4) u>i and s-i axe statistically independent, 

A5) and the perfect knowledge on the noise variance (noise power) o\, is available at the SU. 

B. Preliminaries 

1 ) Energy Detection: In energy detection, the energy of the received signal samples is first computed 
and then is compared to a predetermined threshold. The test procedure of energy detection is given as 

M H t 
i=l H 

where r := [ri,r2, ■ ■ ■ ,tm\, T{r) denotes the test statistic, M denotes the fixed sample size, and 7 e d 
denotes a threshold for energy detection. 

Recall that w,i is a zero mean complex Gaussian RV with variance <7^/2 per dimension. Under Hq, 
the RV 2T(r)/a'f u is a central chi-square RV with 2M degrees of freedom whereas under Hi, the RV 
2T(r)/a^ ] conditioning on |sj| 2 , i = 1, . . . , M, is a noncentral chi-square RV with 2 M degrees of free- 
dom and non-centrality parameter 2\h\ 2 Yli=i \ s i\ 2 l a w- ^ s M increases, 2\h\ 2 ^f=i \ s i\ 2 /°~w approaches 
2M\h\ 2 a 2 /a 2 v , where a 2 denotes the average symbol energy. Let us define SNR m := \h\ 2 a 2 /a^ as 
the minimum detection SNR. Notice that the minimum detection SNR m is a design parameter referring 
to the minimum SNR value by which a detector can achieve the target false-alarm and miss-detection 
probabilities. In practice, it is highly likely that SNR m is different from the exact operating SNR, which 
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is typically difficult to acquire in practice. To distinguish these two different SNRs, we denote by SNR Q 
the operating SNR. 

It follows directly from the central limit theorem (CLT) that as M approaches infinity, the distribution 
of 2T(r)/a^ u converges to a normal distribution given as follows [3] 



2T(r) 



M(2M,AM), undergo, 

(3) 

JV(2M(1 + SNR m ),4M(l + 2SNR m )), under H x . 



Typically, the required sample size M is determined by the target false-alarm and miss-detection 
probabilities, which we denote by a e d and /3 e d respectively. Let M e ™ m be the minimum sample number 
required to achieve the target <5 e d and /3 e d at the detection SNR m level. As shown in [2], we have 



M e 7 n = Ceil(JSNR m )- 2 [Q-^Oed) - Q _1 (l " &d)V2SNR m + lj J = 0(SNR~ 2 ) (4) 

where Ceil(x) denotes the smallest integer not less than x, Q(-) is the complementary cumulative 
distribution function of the standard normal RV, i.e., Q(x) := (27T)" 1 / 2 e~ t2 / 2 dt, and Q~ x (-) denotes 
its inverse function. It is evident from (01) that for energy detection, the number of the required sensing 
samples is inversely proportional to SNR^ when SNR m is sufficiently small [9]. 

As clear from the above description, energy detection has a simple test statistic and has low imple- 
mentation complexity [7]. In addition, known as a form of non-coherent detection, energy detection only 
requires the knowledge on noise power and does not rely on any deterministic knowledge about the 
primary signals Sj. However, one major drawback of energy detection is that, at a low detection SNR 
level, it requires a large amount of sensing time to achieve low detection error probabilities. 

2) A SPRT Based Sensing Scheme: In comparison with a fixed-sample-size detection such as energy 
detection, the SPRT is capable of achieving the same detection performance with a much reduced ASN 
[13]. We next investigate a SPRT based sensing scheme that relies on the amplitude squares of the 
received signal samples [4], [10]. To simplify our description, we now assume that the amplitude squares 
of primary signals, |sj| 2 , i = 1, 2, . . ., are perfectly known at the SU. With this assumption, the spectrum 
sensing problem formulated in £T|) becomes a simple hypothesis testing problem, which is the original 
setting considered by Wald [12]. 

We normalize |rj| 2 as v.- t := 2|rj[ 2 /<7^ for the convenience of derivation. Note that under Hq, Vi is 
an exponential RV with rate parameter 1/2 and under H±, Vi conditional on |sj| 2 is a noncentral chi- 
square RV with two degrees of freedom and non-centrality parameter Aj that can be readily obtained as 
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Aj = 2|/i| 2 |sj| 2 /o"^,. Hence, the probability density function (PDF) of v% under Hq is 

VhM) = \z- Vi/2 (5) 
whereas under H\, the PDF of Vi conditional on Aj is 

p Hl {vi\\) = \e-^ +K)/2 k{V^i) (6) 

where Iq(-) is the zeroth-order modified Bessel function of the first kind. After collecting N samples, 

we can express the accumulative log-likelihood ratio as 

p H (vn\Xn) n n n 

L N (v N \X N ) = log-^ — f f =y^ Zi = - + y^logJ (V^m) (7) 

pv\h (vn) jri tl 

where v N := v 2 , . . . , vn), Zi = log (p Hl (vi\\i)/pH {vi)), and Aat := (Ai, A 2 , . . . , Ajv). The test 
procedure is given as follows: Reject Hq, if Ln(vn\^n) > Accept Hq, if Ljv(«jv| Ajv) < ol; and 
continue sensing, if ol < Lat («jv I ^n) < hi [12], Wald specified a particular choice of the thresholds 
az, and bi for the non-truncated SPRT as follows 

aL = log ^ SP ! t , and 6 L = log 1 _ ^ sprt (8) 

1 — Ctsprt C^sprt 

where Q sprt and /3 sprt denote the target false-alarm and miss-detection probabilities, respectively. For a non- 
truncated SPRT, the Wald's choice on and &l yields true false-alarm and miss-detection probabilities 
that are fairly close to the target ones. 

Let z be a RV that has the same PDF as z\. It has been pointed out in [17] that, in the SPRT, if 
hypotheses Hq and Hi are distinct, then Ejj {z) < < Eh 1 {z), where Eh,(-) denotes the conditional 
expectation under Hi, i = 1, 2. 

As evident from (0, one shortcoming of this SPRT-based sensing scheme is that the test statistic 
contains a modified Bessel function, which may result in high implementation complexity. When the 
perfect knowledge of the instantaneous amplitude squares of the primary signals is not available, the 
PDF under Hi is not completely known, i.e., the alternative hypothesis is composite. Generally speaking, 
two approaches, the Bayesian approach or the generalized likelihood ratio test, can be used to deal with 
such a case. In the Bayesian approach, a prior PDF of the amplitude squares of the primary signals 
is required and multiple summations over all possible amplitudes of the primary signals need to be 
performed, whereas in the generalized likelihood ratio test, a maximum likelihood estimation (MLE) of 
the amplitude squares of the primary signals is needed [18]. Either of these two approaches, however, 
leads to a considerable increase in implementation complexity. 
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III. A Simple Sequential Spectrum Sensing Scheme 
We now present a simple sequential spectrum sensing scheme with the following test statistic 



A 



N 



N 

£ 

i=l 



A) 



(9) 



where A is a predetermined constant. Assuming that the detector needs to make a decision within M 
samples, we propose the following test procedure 



Reject Hq : if K N >b and N < M - 1, or if A M > 7; 
Accept H : if A N < a and N < M — 1, or if A M < 7; 
Continue Sensing : if Ajy £ (a, b) and N < M — 1 

where a, b, and 7 are three predetermined thresholds with a < 0, b > 0, and 7 E (a, 6). 



(10) 

(ID 
(12) 




Fig. 1. The test region of the SSCT. 



In statistical term, the test procedure given in (fT0l- (fT2l is nothing but a truncated sequential test. 
As depicted in Fig. [T] the stopping boundaries of the test region consist of a horizonal line b and a 
horizonal line a, which we simply call the upper- and lower-boundary respectively. Since each term in 
the cumulative sum is a shifted chi-square RV, we simply term the test procedure (fT0l)- ([T2l . the 
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sequential shifted chi-square test. It is evident from (fT0t-(fl2l) that the test statistic depends only on the 
amplitude squares of the received signal samples and the constant A. 

Let r be a RV having the same PDF as |rj| 2 — A, which is the zth incremental term in the test statistic 
©. As we show below, we choose a 2 ^ < A < o 2 w (1 + SNR, n ) to ensure Ejj (r) < < Ej{ 1 (r) similar 
to the SPRT case. In the SSCT, we have E Ho (r) = a 2 w - A and E Hl (r) = o 2 w (1 + SNR D ) - A. Using 
SNR m < SNR G and o 2 w < A < o 2 w (1 + SNR m ), we always have E Ha (r) < < SMR - SNR m < E Hl (r). 
Note that with this choice, the constant A depends on the minimum detection SNR instead of the exact 
operating SNR. 

We normalize the test statistic A at by cr^/2 and rewrite (© as 

N 

Atv = J>< - 2A/a 2 w ) (13) 

i=l 

where A^ := 2Atv/c^ and V{ := 2|rj| 2 /cr 2 y . Let £yv denote the sum of for i = 1,...,N, i.e., 
£at = Yli=i v i an( ^ let ^ denote 2A/cr 2 y . With this notation, we can rewrite A^ as 

An = 6v - NA. (14) 

For notional convenience, we define Ao and £o as zero. Let a, and hi be two parameters defined as 
follows: en = for Nq\ a, = a + iA for i G Np +1 , and 6j = 6 + iA for & G Hq°, where a := 2a/a 2 ] , 
b := 2b/a^, and P denotes the largest integer not greater than —a/A, i.e., P := f loor(— a/A). Using 
£at = X2i=i u *> we rewrite the test procedure (fT0l)-(fT2~l) as 

Reject H : if £ N > b N and N < M - 1, or if £ M > j M ; (15) 

Accept Hq : if ^ < a^ and N < M — 1, or if £jvf < 7m; (16) 

Continue Sensing : if £n G (a.iv,&jv) an d N < M — 1 (17) 

where 7^/ = 7 + MA with 7 := 27/0"^. The corresponding test region is depicted in Fig. [2j where the 
stopping boundaries comprise two slant line segments. We adopt a ssct and /3 ssct to denote the false-alarm 
and miss-detection probabilities of the SSCT, respectively. 

Since the proposed test procedure in (fT0t-(fT2l is not necessarily a SPRT, the Wald's choice on thresh- 
olds, which yields a non-truncated SPRT satisfying specified false-alarm and miss-detection probabilities, 
is no longer applicable. Alternatively and conventionally, the thresholds a, b, 7, the parameter A, and 
a truncated size M are selected beforehand, either purposefully or randomly, and corresponding a ssct 
and /3 ssct are then computed. If the resulted a ssct and /3 ssct do not meet the requirement, the thresholds 
and truncated size are subsequently adjusted. This process continues until a desirable error probability 
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Fig. 2. The test region of the transformed test procedure. 

performance is obtained. In the above process, the key step is to accurately and efficiently evaluate the 
false-alarm and miss-detection probabilities as well as the ASN for prescribed thresholds a, b, 7, the 
parameter A, and a truncated size M, as will be addressed in the following section. 

IV. Evaluations of False-Alarm and Miss-Detection Probabilities 

This section presents the exact false-alarm probability, and a numerical integration algorithm that 
obtains the miss-detection probability in a recursive manner [14]. We start by introducing some preparatory 
tools, including three mutually related integrals that will be frequently used in the evaluation of the false- 
alarm probability. 

A. Preparatory Tools 

The first integral is defined as 

f!g(0 = 1, k = 0, and = f di k d£ fe _! ■■ f 2 rift, k > 1 (18) 
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with Xo = ^ and Xk := [xii • • • > Xk-i,Xk] with < xi < • • • < Xfc- Note that superscript A; and subscript 

(fc) 

Xk are use d to indicate that f Xk '(0 is a fc-fold multiple integral with ordered lower limits specified by 
Xk- Evidently, the integral fxHiO * s a polynomial in £ of degree k. The following lemma shows that 
the exact value of the integral f& (£) can be obtained in a recursive manner (see Appendix |A] for the 
proof). 

Lemma 1: The integral fx^iO is given by 

k— 1 f{k)/£ \k—i 

(fc) 

where the coefficients /• , i = 0, . . . , fc, for fc > 1 can be computed by using the following recurrence 



relation 



fc-i f (fc-l) 

Jk) _ ,(fc-l) . w fc-l ,(fc) _ _Jj / sfc-j 

/i -/i > !eN o > J*, - ~Z^7T37v^ fc -Xi+i) (20) 

i=0 ^ 

(0) (fc) 

and the initial conditions /q = 1. In particular, if xi = X2 = • • • = Xfc = X> /x fc * s gi yen by 

/£ } = ^(e-x) fc . (2D 
Furthermore, the integral /-j? satisfies the following properties 

1) Differential Property: = with Xfc-i = [xi, • • -,Xfc-l] and fc > 2; 

2) Sca/ing Property: = &fj£(€) for f > 0; 

3) SAP P/»/*tfy: /^ } _ 51fc (£ - 5) = /£ } (0- 

It is noteworthy to mention that scaling and shift properties are particularly useful in reducing round-off 
errors when evaluating /jj^ (£). We now introduce the second integral as 

I (0) := 1, and 7 (n) := J ■ ■ ■ J d£ n , n > 1 (22) 

where £ n := • • • with <&<&•••< and fiW = {(6,6, • •.,£„): < & <••• < 

Cn! a i < £i < * ^ ^i }■ I n particular, when n = 1, I*- 1 -* = J 1 d£i = b\ — a\. Let c and d denote two 
non-negative real numbers satisfying < c < d, a^-i < c < b^ and < d. Define the following 
vector 

[»n+l, ■ ■ ■ , On+1, OQ+n+1, • • • , OTV-1, Cj, n fc Nq , 
s .. ' <> ' 



'rn.c 



Q N—Q—n 

[ b n +i,-- - ,b n +i,c ] , n fc N S N 1 Q _ 1 , 

JV-n 

6 n+ l±7V-n, n fc Hf" 2 , 
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where s denotes the integer such that b s < c < b s+ i, Q denotes the integer such that oq < b\ < oq+i, 
and N > 2. Let Aj be an (N — n) x (N — n — i) matrix defined as Ai = [I n "-i-n|0(jv -i~n)xi] T with 
i E ^~ n . We further define the following vectors 



^n,c % = V'n.c • A h i£ Nf n , a£ = [a„ 1+ i, . . . , a n J, n 2 > n x > 



.JV 



yN-n 



(23) 



where i/^" 1 is an (N — i — n) x 1 vector and is an (ni — n<i) x 1 vector. In particular, a™J is defined 
as when n\ = ri2. 

We next show that the exact value of the integral in (1221 can be obtained recursively as follows 
(see Appendix |B] for the proof). 

Lemma 2: The exact value of the integral I^ N ' can be obtained by applying the following recurrence 
relation 

N—2 , \N—n 

/ a « (6jv) " V>2} 2^ (N-n)\ ' 1 



n=0 

Af-2 

rCOrt.^ _ V" AN-n), 



(24) 



/< (^) - e %r ( &jv)/(n) ' 



ri=0 



and the initial condition = 1. 

We now introduce the third integral in the following 



e~^ N d$ 



N 



(25) 



where 9 > 0, N > 1, and := {(&,... ,£at) : < & < . . . < &v; a* < & < h,i E Nf _1 ;c< &v < 
d}. Recalling that c is an arbitrary non-negative number satisfying ajv_i < c < bjy and d is an arbitrary 
number satisfying < c < d and d > a at, we define the following function 



N-n 

n - N e- 9b " +i -E d - l C~ZlJ d > 



N-n 



i=l 



c<h, n E H^~ 2 



/( n) ^ /g^-o^e-fc _ /^:r i} (d) e - 

i=l 

JV-n 

7 (n) J r -iV e -06„ +1 _ J- 0-7^7^^ (d)e- e< 



6d 



, c > bi, n E N 



s-l 




i=l 



c > 61, n E H 



Af-2 



The following lemma shows the exact values of the integral J^P(6) for two particular pairs of (c, d) 
(see Appendix O for the proof). 
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Lemma 3: For any 7^ satisfying < 7a < 6a, the exact values of the integrals «/i^ oo (0) and 



A' 



r(JV-i) 



A-2 



AT 



n=0 



P (A-i), 



-9a„ 



-eb. 



i=l 



(26) 

A-2 
n=0 



We next turn our attention to the evaluation of false-alarm and miss-detection probabilities. 



B. False-Alarm Probability 

Let E N denote the event that A N > b and a < A n < b for n G ttf" 1 with N G and let 

denote the event that Am > 7 and a < A n < 6 for n G Denote by Ph (En) the probability of 

the event En under i^o- Recalling that the test procedure given in ([T0l)-([T2l is equivalent to that given 
in (fT5t-(fl7l). we have 

Ph (ai < & < M G Nf- 1 ; £ N > 6a), N G <~ j 

(28) 

^ P ffo (a, < ii < bi, i G < - 1 ; 6/ > 7m), iV = M. 

Clearly, the false-alarm probability a ssct can be written a ssct = En=i p h„(E N ). In the following 
proposition, we present the exact false-alarm probability a ssct . 
Proposition 1: The false-alarm probability, a ssct , is given by 



Ph (E n ) = { 



M 



(29) 



A=l 



where Ph (En) can be computed as 



A J 



PA 



PA 



PA 



(iV-l)! : 



iV G Nf +1 



A-3 , 



f a N-i (6jv_i)-I{at>3}E' 



6n.+l) 



A-n-1 



A-3 



n=0 



(N-n- 1)1 



-2 n e^P Ho (E n+1 ) 



N G N 



Q+l 
P+2 



n=0 



r M 4 M Ui/2), 



N G K^ 1 
N = M 



'-y M ,oo\ 

where := 2-( Ar - 1 ) e - fejv / 2 . 

Proof: To compute Ph (En), we need to determine the joint PDF of the RVs (£1, ... , £a)- Let 
Py|i/ (fi, • • • , vn\Hq) and P£|# (£i, • • • 1 £aH#o) denote the joint PDFs of the RVs (v±, . . . , vn) and the 
RVs (£1, . . . ,Cn) under Hq, respectively. Recalling that vi is an exponential RV distributed according to 



May 28, 2009 



DRAFT 



IEEE TRANSACTIONS ON SIGNAL PROCESSING (SUBMITTED) 14 

dU, we can write the joint PDF of the RVs (v\, . . . , vn) as Pv\h ( v 1j ■ ■ ■ > vn\Hq) = 2~ N e~ ^*=i Vi l 2 . 
Due to £ N = Ya=\ y i> we nave y i = ^2 = 6 - 6, • • • ) W = £n ~ £n-i, which yields 

Pe\Ho(tu&, ■ ■ ■ , Cn\h ) = pv\ Ho (Si, 6 - 6, • • • , £iv - 6v-i|fli>) = 2- N e-*»P, (30) 

where £o := < £1 < • • • < £at. According to (T28T ) and the definition of T^^, we have 

P Ho (E N ) = P Ho ((t 1 ,b,...,Z N )erW 00 ) = J "... I 2" JV e-^/ 2 ^ Jv . (31) 

b^y , oo 

Generally speaking, because each variable & is lower-bounded by the maximum of ai and and is 
upper-bounded by the minimum of 6j and a direct evaluation is highly complex due to numerous 
possibilities of upper- and lower-limits of (£i,£2, • • • , 6v) [15]. 

Nevertheless, in the case of N E Nf +1 , the parameters a* for i E are all zeros by the definition 

of the parameter P. This implies that ^ is only lower-bounded by for i E 1 , and accordingly 
the upper-bound of £j can be readily identified as 6, t for i E ttf -1 [15]. Using [15, Eqs. (16) and (17)] 
and the fact that is an arithmetic sequence, we obtain Ph (En) as follows 

fbi rb 2 rb N -! roc -b N /2 

P Ho (E N )= dt 2 -- 2- N e-t»/ 2 dt N = (32) 

J & J£l Jb N 2 ~ fJ ! 

with N E ttf +1 . 

We now consider the case of N E #p+2- Since £jv G [&at, oo) and bjy > &jv-i> the upper-limit of 
£at-i is actually b^-i irrespective of £n- Hence, from ((22l . we can write (f3TT > as 

P Ha (E N ) = / 2-^e-W2 d6v ./.../ ^ jy _ i = 2 - {J v-i) e - feN /2 /( iV-i) (33) 

with AT E ^p~2. By applying (g2]|, P Ho (E N ) for iV E N^+J and ^ G ^Q+a* in © can be readily 
obtained. 

We next compute P(Em \Hq). Since 7m E (om,I)m). the upper-limit of £m-i depends on both £jf 
and &Af_i. Thus, the integrals J cZ£a/ and /•••/ are not separable. It is clear from Lemma [3] that 

P(£ M |#o) can be obtained from J^L^) b y setting = 1/2 and N = M in (|26]>. Hence, we have 

P Ho (E M ) = P Ho ((tu ■ ■ ■ ,£m) e Tg^) =/•••/ 2- M e~^/ 2 ^ M = 2" M J^L(l/2). (34) 

T (M) 

Thanks to (|32l ), d33l , and (|34l , we can conclude the proof. ■ 
Remark 1: The problem of evaluating (OTT ) resembles the exact OC evaluation problem in truncated 
sequential life tests involving the exponential distribution [15]. Exact OC for truncated sequential life 
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tests in the exponential case has been solved by Woodall and Kurkjian [15]. However, the Woodall- 
Kurkjian approach [15] is not applicable to evaluate the ASN [16]. More importantly, it cannot be used 
to evaluate d34l ) directly. In the preceding proof, we propose a new approach to derive the exact false-alarm 
probability a ssc t- It is worth mentioning that with slight modifications, our approach is also applicable to 
evaluate the ASN for truncated sequential life tests in the exponential case. 

C. Miss-Detection Probability 

We now turn to the evaluation of the miss-detection probability, /3 ssct . In order to evaluate /? ssct , we 
need to obtain pjj 1 (v{). Since acquiring perfect knowledge on instantaneous Aj or the exact distribution of 
Aj may not be feasible in practice, evaluating PH^Vi) is typically difficult except for the case where the 
primary signals are constant-modulus, i.e., |si| 2 = <r 2 . We next reason that at a relatively low detection 
SNR level, the miss-detection probability obtained by using constant-modulus primary signals can be used 
to well approximate the actual /3 ssct . Our arguments are primarily based on the following two properties 



The first property shows that as N approaches infinity, the distribution of the test statistic £n in the 
SSCT converges to a normal distribution that is independent of a specific choice of Ai, A2, ■ ■ ■ , Aj\r. 
Property 1: The statistical distribution of £/v converges to a normal distribution given by 



as N approaches infinity. 

The property can be readily proved by using the CLT [17]. However, unlike energy detection, this 
property alone is not sufficient to explain that the constant-modulus assumption is valid in approximating 
/3 ssct . This is because each £n for N = I, . . . , M including small values of N, may potentially affect 
the value of /? ssct . To complete our argument, we first present the following definitions. Let £j denote the 
test statistic using the constant-modulus assumption, i.e., |sj| 2 = cr 2 . Define qn '■= bjy for N £ 
and qm '■= 7m for N = M. Let Fjq and Fjy denote the events that £n > Qn, and a, < < bi for 
i G Nf, and £jy > QN and < £j < b{ for i e 'Rf 1 , respectively. Let Ph 1 {Fn) and Ph^Fn) denote the 
probabilities of the events Fn and Fn under H\. Let /3 ssct denote the miss-detection probability obtained 
by assuming constant-modulus signals with average symbol energy a 2 , i.e., |sj| 2 = a 2 . As clear from 
their definitions, we have ssct = YIn=i p hA f n) and /3 ssct = Ejv=i p hA f n)- 

Let A l fi denote the event that Oi <&<&», i € for some integer l N £ , and let A l fi 
denote its counterpart for the constant-modulus case. Let B l £ denote the event that £jv > Qn and 
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&i < ii < bi, i £ an d let B l $ denote its counterpart in the constant-modulus case. We now present 

the second property of the SSCT (see Appendix [D] for the proof). 

Property 2: Let e an arbitrary positive number. If for each N, there exists a positive integer Zjy £ 
such that PhMn) ^ 1 - e/(3M), PhMn) > 1 " e/(3M), and - ^(^)l < e/(3M), 

then 

I Asct — Asct I ^ e 

where depends on the values of N and e. 

Relying on these two properties, we sketch our arguments as follows. To achieve high detection accuracy 
at a low SNR level, the ASN and M are typically quite large. When the sample index N is relatively 
small, it is highly likely that the test statistics £/y an d £n do not cross either of two boundaries. In 
such a situation, there exists some integer I n such that Pjj 1 (A l $ ) and Ph 1 (A l $ ) are fairly close to 1 
whereas Ph 1 (B 1 ] ^) and Ph 1 (B 1 ] ^) are fairly close to 0. Hence, the conditions in Property [2] can be easily 
satisfied. On the other hand, when N is relatively large, one can find a sufficiently large Z/y sucn that 
P Hl (A l £) and P Hl (A l $) are fairly close to one while \P Hl {B l $) - P Hi 0n)\ are sufficiently small due 
to Property Q] guaranteed by the CLT. Collectively, at a low detection SNR level, /3 ssct evaluated under 
the constant-modulus assumption is a close approximation of /3 ssct . Therefore, we will focus on the case 
in which all Aj's are equal to a constant A := 2\h\ 2 a 2 a / = 2SNR m . 

Recall that under Hi, V{ is a non-central chi-square RV, whose PDF involves the zeroth-order modified 
Bessel function of the first kind as given in ©. This makes it infeasible to evaluate /3 ssct by applying the 
computational approach used in the false-alarm probability case. To obtain /3 ssct , we resort to a numerical 
integration algorithm proposed in [14]. 

Defining U{ = V{ — A, we rewrite Ayv in (fT3l) as Ajy = YliLi u i- Clearly, we can write the PDF of Ui 
under Hi as 

PhM) = \e-^ +K+x h Q ^\{u l + A)),u l > -A. (35) 

Recall that M is the maximum number of samples to observe. Denote A-M-k by tk- Let Gk{tk) denote 
the conditional miss-detection probability of the SSCT conditioning on that the first (M — k) samples 
have been observed, the present value t/. = A^_fc, and the test statistic has not crossed either boundary 
in the previous (M — k — 1) samples. If a < t/. < b, an additional sample (the (M — k + l)th sample) is 
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needed. Let u be the next observed value of Uj. The conditional probability Gk(tk\u) can be written as 

if u > b — tk 

1 if u < a - t k (36) 
Gfc-i^fc + u) if a - t k < u < b - t k . 



G k (t k \u) = < 



Using (1361 ). we can compute Gk{tk) as 

ra—tk pb—tk 



/a—tk ro—tk 
p Hl (u)du+ G k -i(t k + u)pH 1 {u)du (37) 

-oo J a—tk 



-oo J a—t k 

for k = 1, • • • , M with the following initial condition 



Go{ta 



(38) 



0, if t > 7 

1, otherwise. 

Employing the above backward recursion process, we can obtain Gm(0), which is equal to the miss- 
detection probability, /3 ssct . 

V. Evaluation of the Average Sample Number 

Roughly speaking, the false-alarm and miss-detection probabilities, and the ASN are three principal 
performance benchmarks for the sequential sensing scheme. In the preceding section, we have only 
concerned ourselves with the error probability performance while in this section, we show how to evaluate 
the ASN. 

Let N s denote the number of samples required to yield a decision. Clearly, N s is a RV in the SSCT, 
and its mean value is the ASN, which can be written as 

E(N S ) = E Ho (N s )P Ha (H ) + E Hl (N s )P Hl (Hi) (39) 

where Eh,{N s ) denotes the ASN conditional on Hi, and Pjj.(Hi) denotes a priori probability of 

hypothesis Hi, for i = 0, 1. 

According to (TTOt — (IT2b. we have 1 < N s < M. Hence, we can express Eh,(N s ) as 

M 

E Hi (N a ) = NP H i {N s = N), i = 0, 1 (40) 

N=l 

where Ph, (N s = N) is the conditional probability that the detector makes a decision at the iVth sample 
under Hi. We now need to determine Pn^Ng = N). Let Cn denote the event that the test statistics 
(£i> • • • > £zv) do not cross either the upper or lower boundary before or at the iVth sample, i.e., 
Cn = £2, • • • , Cat) G b } for N G Hf 1 . For notional convenience, let us define Cq as a universe 
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set. Hence, we have P(Cq) = 1. It implies from (fl3T)-([T7T) that PhA^s = N), i = 0, 1 can be obtained 

as 

P Hi (N s = N) ( = } P Hi (Cjv-i) - (Cat) , iV e - 1 (41) 

P Hi (N s = M) { = } P Ht (C M -x) , N = M (42) 

where the two terms on the right-hand side (RHS) of the equality (a) are the probabilities of the events 
that under Hi, the test statistic does not cross either of two boundaries before or at the (N — l)th sample 
and the iVth sample for N £ _1 , respectively, and the term on the RHS of the equality (b) denotes 
the probability of the event that under Hi, the test statistic does not cross either boundary before or at 
the (M - l)th sample. 

Substituting (0U) and (02]) in (00]), we can rewrite (00]) as 

M-1 M-1 

Eh, (N s )=^N (p Hi (C*_i) - P Ht (Cat) ) + MP Ht (C M -i) = 1 + ^ P Hi (c N ) ■ (43) 

N=l N=l 

According to Lemma [3j OUJ), (0B> and (02]), we have 

ft (^) =2-^,(1/2), iVGKf- 1 . (44) 

Hence, from (|43l , we have 

M-1 

^ (iV s ) = l+E 2 " 7Vj ^L( 1 /2)- (45) 

JV=1 

We next need to evaluate Ph ± {C-n)- Instead of evaluating PhA^n) directly, we compute Ph^n) by 
performing a procedure similar to the one in calculating the miss-detection probability. Note that C C N 
denotes the event that under H\, the test procedure given in (fT5T)-(fT7T) terminates before or at the A^th 
sample, i.e., the test statistic crosses either the upper or lower boundary before or at the iVth sample. 
According to (I38T ). Gk{tk) also depends on 7. With a slight abuse of notation, we rewrite Gk(tk) as 
Gk(tk,j)- Let Vn denote the event that the test statistic crosses the lower boundary before or at the 
iVth sample under H\, and Un denote the event that the test statistic does not cross the upper-boundary 
before or at the iVth sample under H\. It is clear to see Ph x (Vn) = Gn{0, a) and Ph x (Un) = Gn(0, b). 
Obviously, Ph 1 {C c n ) can be written as 

P Hi {Cn) = PhAVn) + (1 - PhAUn)) = G N (0,a) + 1 - G N (0,b) (46) 

where Gjv(t, a) and Gjv(i, &) can be obtained by applying (|37T ) recursively. Hence, we have Ph 1 (Cn) = 
Gat(0,5) — Gjv(0, a). According to (|43l , we have 

M-1 

E Hl (N s ) = 1 + ( G n(0, b) - G N (0, a)). (47) 

N=l 
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In the following proposition, we present the ASN of the SSCT. 

Proposition 2: The ASN of the SSCT can be obtained as 

Af-l M-l 

E(N S ) = P Ho (H ){l+Y J 2-"jW>„(l/2)) +P Hl (H 1 )(l+ £(Gfr(0,5)-G w (0,a))). (48) 

N=l N=l 

Proof: The proof follows immediately from d39l ), d45l ), and (|47T ). ■ 

VI. Simulations 

In this section, we provide several examples to illustrate the effectiveness of the SSCT. In all simulation 
examples, we select A to be 2 + SNR m , which ensures Eh (An) = —Eh 1 (An). In the first three test 
examples, the truncated sample size M is selected to be the minimum sample number M e ™ m required by 
energy detection to achieve specified false-alarm and miss-detection probabilities, and a is always chosen 
to be —b. All the test examples assume that \h\ is equal to one, and adopt constant-modulus quadrature 
phase shift-keying (QPSK) signals except for Test Example 2, in which the modulation formats of the 
primary signals are explicitly stated. Following the conventional terminology in sequential detection, we 
define the efficiency of the SSCT as r) := 1 - ASN/M e ™ in . 

Test Example 1 (The SSCT Versus Energy Detection): Table U compares the SSCT with energy detection 
in terms of false-alarm and miss-detection probabilities and the ASN for different SNR m . The truncated 
sizes corresponding to SNR m = 0, —5, —10, and —15 dB, are selected to be the minimum sample sizes 
required by energy detection to achieve target (a ssct , /3 ssct ) = (0.01,0.01), (0.05,0.05), (0.1,0.1), and 
(0.15,0.15), respectively. The parameters b, A and 7 are given in the table. It is shown in Table U that 
compared with energy detection, the SSCT can achieve about 29% ~ 35% savings in the average sensing 
time while maintaining a comparable detection performance. It can be also observed from the table that, as 
indicated by using an abbreviation, Numerical, in the parenthesis, the false-alarm probabilities computed 
from the exact result (l30l and the miss-detection probabilities obtained by the numerical integration 
algorithm match well with those obtained by Monte Carlo simulations. 

Test Example 2 (Detection Performance Without Knowing the Modulation Format of the Primary 
Signals): In this example, the primary signals are modulated by using a square 64-quadrature amplitude 
modulation (QAM) with a 2 s = W sm ™/ W al . Table HI compares with the miss-detection probabilities and 
the ASN between the SSCT and energy detection for SNR m = 0, —5, —10, and —15 dB. The parameters 
a, b, 7, A, and M in the SSCT, are determined by using constant-modulus QPSK signals with the 
same average symbol energy a 2 = 10 SNRm//10 cr^, and these design parameters are used in the SSCT to 
detect the 64-QAM primary signals. That is, the SSCT does not have the knowledge of the modulation 
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format of the primary signals. When evaluating the miss-detection probability and the ASN of the SSCT 
in the 64-QAM case, we use equiprobable prior distributions of Sj to obtain pfj 1 (ui). As can be seen 
from the table, the miss-detection probabilities obtained in the 64-QAM case are fairly close to the ones 
obtained in the QPSK case except for the case of SNR m = dB, where M = 40. This is because M is 
not large enough to neglect errors caused by using the CLT approximation. However, energy detection 
and the SSCT suffer from a similar amount of approximation error when evaluating the miss-detection 
probability. 

Test Example 3 (Mismatch between SNR Q and SNR m ): Table UII1 lists the false-alarm and miss-detection 
probabilities and ASN when SNR D = -12, -13, -14 dB, and SNR m = -15 dB. The parameters for the 
SSCT and energy detection are determined to ensure that a target pair (a ssct , (3 ssct ) satisfies (0.15,0.15) 
at SNR m = —15 dB. It can be seen from the table that, as SNR G increases, the miss-detection probability 
decreases while the false-alarm probability remains unchanged. All obtained (a ssct , /3 ssct ) pairs satisfy the 
target false-alarm and miss-detection probability requirement. Interestingly, the efficiency of the SSCT 
increases to 46% from 29% even though the miss-detection probabilities of the SSCT are larger than 
those in energy detection. It implies from this example that the SSCT offers more flexibility in striking 
the tradeoff between sensing time and detection performance than energy detection. 

Test Example 4 (Impacts of Truncated Size M on the Efficiency of the SSCT): Let T p denote the 
probability of the event that the SSCT ends at the Mth sample (truncated at the Mth sample). Table 
JV] lists T p , the ASN, and the efficiency of the SSCT for various selected combinations of M, a, b, and 
7 at SNR m = —5 dB to achieve target (a ssct ,/? S sct) = (0.055,0.046). The results shown in the table are 
obtained by using Monte Carlo simulation. To achieve roughly the same false-alarm and miss-detection 
probabilities, the sample size for energy detection is chosen to be 140. As can be seen from this table, the 
efficiency of the SSCT increases as the truncated size M increases but the pace of the improvement is 
diminishing. Table UVl also lists the efficiency of the non-truncated SPRT based sensing scheme presented 
in Section III-B2I It is clear from the table that as the truncated size M increases, the efficiency of the 
SSCT comes fairly close to the one achieved by the non-truncated SPRT based sensing scheme. 

VII. Conclusion 

In cognitive radio networks, stringent requirements on the secondary and opportunistic access to 
licensed spectrum, necessitate the need to develop a spectrum sensing scheme that is able to quickly 
detect weak primary signals with high accuracy in a non-coherent fashion. Motivated by this, we have 
proposed a sequential sensing scheme that possesses several desirable features suitable for cognitive radio 
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networks. To efficiently and accurately obtain major performance benchmarks of our sensing scheme, 
we have derived an exact result for the false-alarm probability and have applied a numerical integration 
algorithm to compute the miss-detection probability and the ASN. 

There are several potential extensions of this work that deserve further exploration. First, our approach 
to determining design parameters such as thresholds and the truncated size follows the original Wald's 
approach in the sense that the cost of observations as well as the cost of the false-alarm and miss-detection 
events have not been considered. A Bayesian formulation of the SSCT can be an interesting extension. 
Second, this work assumes perfect knowledge on noise power at the SU, which perhaps is difficult to 
acquire in practice. The effect of noise power uncertainty on the SSCT is worth investigating. Third, 
another extension of this work is to study sensing-throughput tradeoffs for the SSCT. 

Appendix A 
Proof of Lemma [T] 

We prove the lemma by induction. It is obvious from ( fT8l) that ( fT9l holds for k = 1. Now suppose 
that ( fT9l and (l20b hold for the case of k — 1. By definition and the induction assumption for k — 1, we 
have 



r? /k-2 ,(fc-l) \ 



i 



k-2 Ak-l) 
J i 



= E 7hv&-i " + fti ■ « " Xk) 

i=0 {K-V- Xk 

fc-l -(fc-l) fc-2 Ak-l) 

- £ fihji« - - £ (b^ - m 

i=o y ' i=0 K 1 

Clearly, comparing ( fT9l with d49l , we can readily conclude the recurrence relation given in (l20b . In 

particular, when \i = X2 = • • • = Xfc> au coefficients except are zeros and hence (|2TT) follows 

immediately. 

Since the differential property can be proved in a straightforward manner, we omit the proof. We next 
prove the scaling property by induction. When k = 1, we have 



'iXi 

Hence, the scaling property holds for A; = 1. We now suppose that the property holds for k = n — 1. 
Applying the induction assumption and a substitution = u, we can rewrite the integral fjjr (^) as 
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Hence, the scaling property holds for k = n. This concludes the proof. The shift property can be proved 
in a similar manner and thus the proof is omitted. ■ 

Appendix B 
Proof of Lemma [2] 

Recall that £jv is a sum of Vi, i.e., £jv = Yli=i v i with Vi = 2\n\ 2 /a^,. Hence, (£1, ■ ■ ■ , Cat) is 
a non-decreasing sequence, i.e., £o < £1 < £2 < ■ ■ ■ < £a- Fig. |3(a)| plots the parameter ^tv versus 
the sample index iV. It is clear from its definition that the region contains all possible sequences 
(£i> • • • , £at) (simply called paths hereafter) satisfying £0 < £1 < £2 < • • • < £a and < a,, < & < 6j. 
Hence, the zth component of each path (£i,£2, • • • 5 £a) in is lower-bounded by the maximum of 
and a^, and is upper-bounded by the minimum of £j + i and The direct computation of I^ N > 
is highly complex [15] due to numerous possibilities for lower- and upper-limits in the integral I^ N \ 
Considering the fact that the integral can be readily computed if either the lower- or upper-limit is 
a constant, we express into an equivalent set, over which the integration can be readily computed 
in a recursive fashion, thereby obviating the need to exhaustively enumerate these possibilities. 

Let (pi, i € N^, denote a sequence of real numbers with < cj)\ < 4>2 < . • • < 4>n- Let us first define 
the following set, 

^ = {(£ n+ i,£ n+ 2, • • • , £a) : 

< £ n +i < £n+2 . . . < &v; fa < £i < fi+i, i G i^i 1 ; 0a < £a < 6a} (50) 

where 0„ := [<^ n +i, . . . , 4>n] is an (N — n) -dimensional real vector with the ith entry of the vector cf>^ 
being the lower bound of £ n+ j for i £ ^~ n . We next define the following non-overlapping subsets of 

U (N) 



»W := {(&, . . . ,£a) : (6, • • • ,£at) G n^ } , ai <Zi<bi,ie N?; 6 n+1 < £ n+1 < £„+ 2 } (51) 
where n G H^~ 2 - 

As can be seen from Figs. |3(b)| and |3(c)[ 11^/ contains all possible paths, (£i,£2, . . . , £a), which 
are lower-bounded by (01, 02, . . . , a at) and upper-bounded by (£2, £3, ■ ■ ■ , £a, &a)> whereas for 
n G N^ -2 contains all possible paths (£1, £2, • • • 5 £at) having the property that the first n variables lie in 
the set i.e., (&,...,&,)€ ft (n) , and the (n + l)th variable £ n +i excesses the upper slant line, i.e., 
£n+i > 6n+i- Again, it is clear from its definition that the set £l( N > is equal to the difference between 
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(c) Possible Paths from E^'. (d) Possible Paths from Q (N) (N > Q). 



Fig. 3. An illustration for Proof of Lemma [2] 



the set IL^J and the union of E^ for n G H^~ 2 , i.e., = U^J n ( U^~ 2 E { n N) ) c . Thus, it follows 
from d22]) that 

N— 2 

i {N) = /•••/'^v (52) 



rr 



n=0 



We now evaluate two terms on the right-hand side (RHS) of (I52I ). It is clear from (1181 ) and (1501 ) that 
the first term on the RHS of (1521 ) is nothing but / ^ (&zv)- We next take a close look at the second term. 
The evaluation of the second term is categorized into the following two cases: 

• Case 1: N < Q: In this case, we have b n+ \ > b\ > a n for any n € {1, . . . , N}. It implies from 

(ED that we can express E^ as E^ = (n) x [b n+1 ,£ n+2 ] x ••• x [6 n+1 ,^] x [b n+1 ,b N ] for 
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n = 0, . . . , N — 2 where := 0. According to (l50l ). we have 

= n< n) x <X n , nG <" 2 - (53) 

Since £ n+ i > b n+ \ > b n > £ n , the integral over S7( n ) and that over n^~^_ are separable. Hence, 
relying on (|24l) . (f5Tb and (|53T ). we have 

/ "' I d ^ N = I "' / d ^+ i ---^ x / ••• / ■• (M^ (n) - (54) 

• Case 2: N > Q + 1: The proof in this case follows the same line of argument as that in the previous 
case. The key difference is that because N > Q + 1, some a n may be larger than b n+ \, as depicted 
in Fig. |3(d)| To be specific, from the definition of Q, we have 

a-Q < b\ < aq+i, . . . , aq +n < b n+ i < a,Q +n+ i, . . . ,ajv < &at-q+i < ojv+i- 

1) For n G H^ - ^" 1 , we have 

~n = ^ (n) x[6 n+ i,^ n+2 ] x ••• x [6 n+ i,^ n+ Q+i]x[aQ +n+ i,^Q +n+2 ] x ••• x [ajv-i,^] x [a N ,b N ]. 

V v ' V v ' 

Q N-Q—n 

Equivalently, = ft( n ) x n^~ n) with tj)% = [6 n+1 , . . . , 6 n+1 , flQ+n+i, • • • , ajv]- 

v„, ajv > ™ v „ / > v / 

Q N-Q-n 

2) For n G ^jv-Q' we have < b n+ \. Due to o« < aw for i G N^ -1 , we have a-i < b n+ \ for all 
i G K^ -1 . Since any (£ n +i 5 • • • , Cat) 6 n^t _n ^ belongs to a Cartesian product of (N — n) intervals 
(a hyper-rectangle) [b n+1 ,£ n+2 ] x ... x [&„ + i,£at] X [6„+i,&jv] having the same lower limit b n +i, 
this case is the same as Case 1. Equivalently, = fifa) x LT^ n ^ with tpn,a N = b n +ilN~n- 
Summarizing the preceding results for Case 2, we have 

/ ■" I d$N= I ■" / dfn + i ""^ x / ••• / ^•••^ =/ s~r )(67v)/(n) (55) 

where ^ = [6 n+ i, . . . , 6 n+1 , OQ +n+ i, . . . , ajv] for n G H^~ Q_1 and = 6 n+ ilAr_„ for 

V v ' > „ ' 

Q N-Q-n 

n G H^" 2 



The proof follows immediately from (T52l . (I54I ). and (I551 ) and the fact that the first term on the RHS 
of ([52]) is equal to fffi(b N ). ■ 
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Appendix C 
Proof of Lemma [3] 

Though the idea of the proof can be extended to a general case of a^-i < c and a at < d, we will 
consider the following two cases: Case 1: c := 7at < and d := oo, and Case 2: c := a at and d := b^, 
which correspond to (l26l) and (l27l) respectively. Define the following sets 

©S } : = {(6, ... :&<&<■■•< 6v; Oi < & < fc+i, * e <~\ c < £v < d}, (56) 
*S' n) : = • • • M : ■■■^N)e e^JOj < & < M € ^;6 n+1 < < £n +2 }, (57) 

where are non-overlapping subsets of f° r n £ ^o^ -2 - The integral over Q^ N J can be readily 

computed as 

[■■■ [ e- e t»d$ N = /"V^d&v f d^N-i--- t\i 

N 

= J2 d ~V ( J-Hc)e- ec - f$Z!\d)e- dd ) (58) 

* — * «*o "o 

i=l 

where the equality in (|58l ) is obtained by using integration by parts repeatedly and the differential property 

4f£ ) (o/de = /£- 1) (o. 

• Case 1: oat < 7at < ftjv and d = oo. Similarly to the argument used in Lemma |2j we have 

Tgjoo = <!oo n ( Ul" 2 <; n 2) C and thus we have 

N— 2 

<L(*) = /•••/ - E / ■ / e "^^iv (59) 

Substituting c = 7at and d = oo in (|58l ) and using the fact that e _£,d is zero for 9 > 0, we have 

AT 



[■■■ [e-^dt N = Y J 0- l f i a p\j N )e- t 
i=i 



1) For 7at < 6i, we have 6i > since 77V > otv- Since £ n < fo n < b n+ \, we have ^^'^ = X 
Il^ N ^ and the integrations over ft( n ) and nr ? J are separable. Thus, applying integration 
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by parts and the fact fz? k (xfc) = 0. we obtain 



e~ e ^di n+l ---di N x 



n 



fiC* 



j(n)gn-N e -eb n + 1 _ 



(60) 



2) For 7tv > &i, we have s > 1 since 6 S < 7/v < b s +i- Similarly to the argument used in Case 2 of 
the proof of Lemma H we have $W = 0( n ) x n%" n) , 



A? 



^A-n JN-n-i) 
t(N-n-i), 



n G N 







/( „) E iV-n ^ 7 C;v :r -) (6n+l)e -^ +1) „ e ^A-2 



(61) 



This concludes the proof for Case 1. 

Case 2: c = a/v and b = bjy. Similarly to the method used in Case 1, we have T^ J 
©£L n ( <" 2 <W and thus we can express ^.(fl) as 



(AT) 



N-2 



a N ,b N a N ,b N 



'd£ 



A- 



(62) 



The rest of the proof is analogous to that in Case 1. The key difference is that in Case 2, the term 
e -6d j s no i on g er Z ero. From d58l ), the first term on the RHS of (|62l can be readily obtained as 

N 



(63) 



G 



i=l 



Since 6tv_q < ajv < b^-Q+i, we have s = N — Q in this case. Similarly to the method used to 
derive (l60l and (|6TT >. we have 



-66 



A 



*(»'"' 

a N ,b N 



N-n 



i=i 



>AT-2 



N-n 



j(n) ^ 0-* 



i=l 



(N-n-i) 



N-n 



i=l 

The proof for Case 2 follows clearly from (l62l . (l63l . and (l64l ). 



®N < G ^0 
OA? > 6i,n € 



(64) 
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TABLE I 

SCCT Versus Energy Detection 



SNR m (dB) 





-5 


-10 


-15 


7 


-8.5 


-5.69 


-4 


-1.897 


b 


27 


35.32 


69.30 


158.47 


A 


3 


2.316 


2.100 


2.032 


Qssct (Monte Carlo) 


0.011 


0.055 


0.103 


0.153 


a ssct (Numerical) 


0.011 


0.055 


0.103 


0.153 


a e d (Energy Detect.) 


0.011 


0.055 


0.101 


0.150 


/3 SS ct (Monte Carlo) 


0.008 


0.046 


0.099 


0.154 


/3 ssct (Numerical) 


0.008 


0.047 


0.100 


0.156 


/3 c(i (Energy Detect.) 


0.008 


0.046 


0.096 


0.149 


ASN (Monte Carlo) 


26 


95 


509 


3154 


ASN (Numerical) 


26 


96 


515 


3185 


M (Energy Detect.) 


40 


140 


730 


4450 


Efficiency 77 


35% 


32% 


30% 


29% 



Appendix D 
Proof of Property [2] 

Note that F N = A l $ nB l fi and F N = A l $ DB 1 ^. Applying the inclusion-exclusion identity [19, p. 80], 
we have P Hl (F N ) = P Hl {A l £ ) + P Hl ) - P Hl (Aft u B N ) and p H, ( p n) = P Hl (Aft ) + p ft {&$ ) " 
Ph 1 (A 1 ^ U B°fi). Thus, by using the triangle inequality, we have 

\Pm(F N ) - P Hi (F n )\<\P Hi (A 1 n) ~ p hMn n )\ + \ p Hl (B L N N ) - PhAB L n)\ 

+ | P Hl (Afr U B^f) - P Hl {A L N N U B L N » ) | . (65) 

Since 1 - e/(3M) < P Hl {A L N N ) < 1 and 1 - e/(3M) < P Hl {A L N N ) < 1, we have 1 - e/(3M) < 
PhMn n ) < p hMn n U B% n ) < 1 and 1 - e/(3M) < PhA^n) < p hMn n U B^ 1 ) < 1. 
It can be readily inferred from the above inequalities that \P Hl (A^ N ) - P Hi (-^n n )\ - e/(3M) and 
\P Hl (Ax N UB^ n )-Ph 1 (An n U ®N N )\ ^ e/(3M). This, along with the inequality (gSJ) and the assumption 
\P Hi {I3n n ) ~ p hABn n )\ < e/(3M), implies that \P Hi ( p n) - P Hi ( p n)\ < e/M. Hence, we have 

M 

lAsct - Asctl < \ P h,(Fn) ~ p hAFn)\ < e. 
i=i 
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TABLE II 

Detection Performance Without Knowing Modulation Types of the Primary Signals 



SNR m (dB) 





-5 


-10 


-15 


Asset (QPSK,Monte Carlo) 


0.008 


0.046 


0.099 


0.154 


Asset (QPSK,Numerical) 


0.008 


0.047 


0.100 


0.156 


/3ed (QPSK,Energy Detect.) 


0.008 


0.046 


0.096 


0.149 


Asset (64-QAM,Monte Carlo) 


0.012 


0.048 


0.099 


0.154 


Asset (64-QAM,Numerical) 


0.012 


0.050 


0.103 


0.157 


Asset (64-QAM,Energy Detect.) 


0.012 


0.047 


0.096 


0.149 


ASN (QPSK,Monte Carlo) 


26 


95 


509 


3154 


ASN (QPSK,Numerical) 


26 


96 


515 


3185 


ASN (64-QAM,Monte Carlo) 


26 


95 


509 


3154 


ASN (64-QAM,Numerical) 


26 


96 


514 


3190 



TABLE III 

Mismatch between SNR m and SNR q (SNR m = -15 dB) 



SNR (dB) 


-12 


-13 


-14 


-15 


Asa (Monte Carlo) 


0.0018 


0.0153 


0.0629 


0.154 


Asct (Numerical) 


0.0017 


0.0151 


0.0628 


0.156 


Aed (Energy Detect.) 


0.0012 


0.0131 


0.0584 


0.149 


ASN (Monte Carlo) 


2425 


2686 


2948 


3154 


ASN (Numerical) 


2499 


2769 


3035 


3185 


M (Energy Detect.) 


4450 


4450 


4450 


4450 


Efficiency (ry) 


46% 


40% 


34% 


29% 



TABLE IV 

Impacts of Truncation Size M (SNR = -5 dB, Monte Carlo Simulation) 



M 


a 


b 


7 


ASN 


T p 


V 


M = 140 


-35.32 


35.32 


-5.69 


95.4 


26.8% 


32% 


M = 160 


-28.95 


23.16 


-5.50 


76.7 


9.2% 


45% 


M = 180 


-27.33 


21.54 


-6.00 


73.1 


5.1% 


48% 


M = 200 


-26.40 


20.85 


-6.32 


71.2 


3.0% 


49% 


M = 500 


-25.48 


19.69 


-6.32 


68.8 


0.005% 


51% 


M = 1000 


-25.42 


19.63 


-6.32 


68.6 





51% 


SPRT (non-truncated) 








67.9 





52% 
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